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In Brief 
Molecular basis for recognition of exogenous agonists by the human bombesin 
receptor subtype-3 provides avenues for design of new modulators for diabetes and 


obesity. 


Highlights 
e Structures of BRS3-G, complexes, both in the absence (apo-form) and 
presence of the ligands BA1 and MK-5046, were determined. 
e The important residues in the ligand binding pocket for bombesin receptor 
subtype-3 activation were described. 
e Тһе basis for ligand selectivity of bombesin receptors was depicted in detail. 
e Resolution of the binding pockets provides path for rational drug design against 


bombesin receptor subtype-3 for the treatment of diabetes and obesity. 
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SUMMARY 


Bombesin receptor subtype-3 (BRS3) is an important orphan G protein-coupled 
receptor that regulates energy homeostasis and insulin secretion. As a member of the 
bombesin receptor (ВпК) family, which includes neuromedin В receptor (NMBR) and 
gastrin-releasing peptide receptor (GRPR), the lack of known endogenous ligands and 
high-resolution structure has impeded understanding of BRS3 signaling and function. 
Here, we present cryogenic electron microscopy (cryo-EM) structures of BRS3 in 
complex with heterotrimeric G, protein in three states: apo, bound to the pan-BnR 
agonist, BA1, and bound to the synthetic BRS3-specific agonist MK-5046. These 
structures reveal the architecture of the orthosteric ligand pocket underpinning 
molecular recognition. Comparisons with BnR members provide insights into the 
structural basis for BRS3’s selectivity and low affinity for bombesin peptides. 
Examination of conserved micro-switches suggests a shared activation mechanism 
among BnRs. Together our results enable deeper exploration of BRS3’s ligand 


selectivity, signaling, and therapeutic targeting for diabetes and obesity. 


Introduction 


G protein-coupled receptors (GPCRs), pivotal membrane proteins, orchestrate cellular 
responses to a variety of extracellular stimuli, including peptide ligands. Among these, 
bombesin, a 14-amino acid peptide originally isolated from the European xenopus 
laevis (Bombina bombina)!?, stands out. In mammals, bombesin analogs such as 
neuromedin B (NMB) and gastrin-releasing peptide (GRP) interact with specific 
receptors in the bombesin receptor (BnR) family?. These receptors are integral to 
physiological processes in the gastrointestinal tract and central nervous system, 
influencing aspects like food intake^, fear memory consolidation?, and itch sensation". 
Structural studies of NMBR and GRPR with various ligands have shed light on the 
interactions within the BnR family, enhancing our understanding of their structure- 
function relationships and aiding in the development of novel therapeutics for BnR- 


related conditions, including cancers and pruritus*?. 
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Bombesin receptor subtype-3 (BRS3), also a member of the BnR family, is an orphan 
GPCR primarily found in the brain. Unlike its counterparts NUBR and GRPR, BRS3 
does not exhibit high affinity for the known endogenous bombesin family 
peptides*"°"(Table S1), and its natural ligand remains elusive. Despite this, specific 
ligands have been identified for BRS3. Among them, the synthetic bombesin analog 
[D-Phe$, g-Ala", Рһе!3, NIe'^]Bn(6—14) (where D and В indicate Dextro- and В- 
conformations, respectively), referred to as BA1, exhibits high affinity, binding to all 
three bombesin family receptors. As a pan-BnR agonist, BA1 is frequently utilized in 
research as а BRS3 agonist/?'?. Nonpeptide agonists like MK-5046 have shown 
selective BRS3 agonism, enhancing our ability to investigate BRS3's functions 1415, 
Additional agonists, such as DSO-5a and oridonin, have displayed potent and selective 
BRS3 activation, suggesting their potential in treating metabolic disorders 1977, 
Particularly, in comparison with MK-5046, DSO-5a activates peripheral BRS3 to control 
blood glucose homeostasis while circumventing central effects ^. Additionally, oridonin, 
derived from Isodon rubescens and commonly employed in Traditional Chinese 
Medicine", exhibits pharmacological activities encompassing hypoglycemic and anti- 
inflammatory effects'?, as well as anti-cancer effects?) As a novel BRS3 ligand, 
oridonin holds promise as a lead compound for treating metabolic disorders, including 


type 2 diabetes and obesity. 


Utilizing these synthetic ligands and animal models has revealed BRS3's critical roles 
in regulating energy homeostasis, glucose metabolism, and blood pressure. Its 
significant function in obesity was initially дел ед”, followed by findings underscoring 
its necessity for energy balance and insulin secretion?^?^, BRS3 neurons in the 
hypothalamus can be activated by cold exposure, refeeding, and certain agonists. 
Notably, activation of BRS3 neurons in the dorsomedial hypothalamus increased 
sympathetic nervous tone, elevating energy expenditure, body temperature, heart rate, 
and blood pressure via brown adipose tissue, heart, and arteries, respectively, without 


affecting food intake or physical activity. Conversely, the activation of BRS3 neurons 
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in the paraventricular nucleus of the hypothalamus reduced food intake, without 
affecting energy expenditure, body temperature, or physical activity? (Figure1A). 
Additionally, agonists of BRS3 can stimulate insulin secretion from pancreatic 8 cells 
(Figure 1A), making BRS3 a promising target for metabolic diseases such as obesity 
and diabetes? ?. However, the lack of structural data for BRS3, especially in different 
ligand-bound states, has hindered our understanding its molecular mechanisms and 


the development of targeted treatments. 


In this study, we present cryo-EM structures of BRS3 in active conformations, both in 
the absence and presence of ligands (BA1 and MK-5046). Together with mutational 
studies, these structures elucidate the mechanisms of ligand binding and receptor 
activation, offering crucial insights into this orphan receptor. This information is 
invaluable for structure-based drug design targeting BRS3, potentially revolutionizing 


treatment approaches for obesity, diabetes, hypertension, and related diseases. 


Results 


Structure determination of BRS3-G, complexes 

For the cryo-EM analysis of BRS3-G4 complexes, we engineered the human BRS3 by 
attaching a thermostabilized apocytochrome b*°®RIL (BRIL) to its N-terminus. This was 
done to enhance the stability and homogeneity of the complex, crucial for effective 
cryo-EM studies?’ (Figure S2). Additionally, we utilized the NanoBiT tethering strategy: 
fusing the large part of NanoBiT (LgBiT) to the C-terminus of BRS3 and a 13-amino 
acid NanoBiT peptide (HiBiT) to the C-terminus of the GB subunit??. We also employed 
an engineered G, protein chimera, similar to the one used in prior structural studies of 
NMBR and GRPR®. Here, "Са refers to this specific engineered chimera. We 
assembled the BRS3-G, complexes for cryo-EM by co-expressing and purifying BRS3- 
LgBiT, Са Gfy-HiBiT, antibody scFv16, and nanobody-35 (Nb35)?99 These 
components were processed both in the absence (apo-form) and presence of the 
ligands BA1 and MK-5046. This approach facilitated the efficient formation and 


enhanced stability of the apo and ligand-bound BRS3-G, complexes. 
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The single-particle cryo-EM analysis of the collected datasets yielded density maps at 
nominal global resolutions of 3.5, 2.9, and 3.0 A for the apo BRS3-G,-scFv16, BA1- 
bound BRS3-G, and MK-5046-bound BRS3-G,-scFv16-Nb35 complexes, 
respectively (Figures 1B-1G, S3-5; Table S2). These maps clearly delineated the 
positions of the receptor, G-protein trimer, and scFv16/Nb35 in each complex, with 
most amino acid side chains being resolved (Figure S6). In the apo BRS3-G,-scFv16 
complex, the orthosteric binding pocket (OBP) did not exhibit any additional densities. 
However, for the BA1- and MK-5046-bound complexes, densities for BA1 and MK- 
5046 were well-resolved, allowing for the construction of atomic models with a high 
degree of confidence. All three complexes displayed a similar overall structural 
arrangement, with root mean square deviation (RMSD) values for Ca atoms ranging 


between 0.3-0.6 À, indicating minor variations among them (Figures S7A-S7C). 


Recognition of BA1 by BRS3 

In our study, the quality of the cryo-EM maps allowed for the precise modeling of all 
nine amino acid residues of ВАТ. This finding aligns with our previous structures of 
NMB-NMBR and GRP(14-27)-GRPR, where only the C-terminal 9 residues of the 
longer bombesin agonists (33-residue peptide NMB30 апа 14-residue peptide 
GRP(14-27), respectively) were observed, despite the utilization of the 33-residue 
peptide NMB30 and the 14-residue peptide GRP(14—27) during complex preparation, 
respectively (Figures 2A and 2B; Table S3). This suggests that a minimal 9-residue 
segment might be sufficient for full agonist potency. In the BRS3 structure, BA1 is 
positioned within the central OBP located in the upper transmembrane domain (TMD), 
interacting with all extracellular loops and transmembrane segments except TM1. The 
peptide's C-terminus penetrates deeply into the helical bundle, while the N-terminus of 
BA1 extends outward towards the extracellular ligand binding cavity. BA1 assumes a 
dumbbell-shaped conformation, reminiscent of NMB and GRP(14-27) in their 


respective complexes, underlining a conserved recognition pattern among bombesin 
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receptor family agonists (Figure 2B). For clarity, all positions of BA1 are numerically 


designated from the amino terminus of NMB or NMC (Figure 2A). 


The nine residues of BA1 (residues 2-10) can be categorized into three motifs: the C- 
terminal НЕХ(8-10)84 motif (Where the X is Ме), the middle АУВА(5-7)8 motif, and the 
N-terminal PFQW(2-4)®*' motif (Figure 2A). Тһе HFX(8-10)®4' motif adopts ап 
extended structure at the bottom of the ОВР, with the side chain of X10®4' penetrating 
the deepest. The side chains of Х108М and F9?^! are encased іп a hydrophobic cavity 
comprising 16253 1128333 V131??6, (221542 | 225546 W284°48, A319742, and 
320743 (Figure 2C). These hydrophobic interactions are crucial for BRS3 activation, 
as demonstrated by significantly reduced potency when 196253 1128333, апа М131336 
are substituted with alanine (Figures 2F and S8; Table 54). In addition, F9P^' engages 
in stabilizing п-п packing interactions with the surrounding Е225548, H28895* and 
Y29185*(Figure 2D). The H8?^' residue forms hydrogen bonds with 8124329 and 
5205502 differing from the sandwiching of the corresponding H8 residues in ММВ and 
GRP(14-27) by P??? and PEC! in the NMB-NMBR and GRP(14-27)-GRPR structures 
(Figures 2D and S7F). The alanine mutations of 5124329 ог 5205-X2 had a minimal 
effect on BRS3 activation by ВАТ, yet, the proline substitution of the 5124329 or 
5205-X2 substantially lowed the BRS3 activity coordinated by ВАТ (Figures 2D and 
S8; Table S4). This highlights the pivotal roles of 5124329 and S205®°' in peptide BA1 
recognition by BRS3. H8? also forms hydrogen bonds with К127332 and 
E182^9*(Figure 2D). Furthermore, the К127332 also forms a hydrogen bond with the 
backbone carbonyl of the F95, and the R316’°° forms a hydrogen bond with the 
backbone carbonyl of the H8? (Figure 2D). The carboxyl group of the Х1ОВАТ forms two 
hydrogen bonds with the side chains of the К127332 and КЗ16739 (Figure 2D). Mutation 
of R127?*?, E1824, С221542 Ң288652 Y29155 and R316"?? in BRS3 results in a 
reduced BA1-induced signaling response (Figures 2F, 2G and S8; Table S4), 


confirming the essential roles of these residues in BA1 binding and receptor activation. 
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Тһе AVPA(5-7)?^! motif primarily engages in hydrophobic interactions with F191E°! and 
Е200ЕС2 from the ECL2, and P304*^? from ECL3 (Figure 2D). Indeed, alanine 
mutation of the Е200ЕС22 resulted in decrease in Са signaling (Figures 2G and S8; Table 
S4). In addition to the hydrophobic interactions, the backbone carbonyl group of PA7P^! 
formed a hydrogen bond with the backbone amide group of S205*€"?, further stabilizing 
ВАЛ binding. 


Following the middle motif, the PFQW(2-4)9^! motif extends along the upper pocket, 
engaging in T-T interactions with Y 1082955, F200*€72, Ң308731 and F309"?? from BRS3 
(Figure 2E). Alanine substitutions at these sites diminish receptor activation compared 
to the wild-type BRS3 (Figure 2G; Table S4). Апоп-п interactions between °F2°4' апа 
D193**2/D303**P? also contribute to BA1’s binding stability, with alanine mutation of 
D303*€Ó significantly reducing Са signaling (Figure 2G; Table S4). Furthermore, Q3B^1 
forms extensive hydrogen bonds with the 8:53 backbone, and its carbonyl group 
interacts with H1072% (Figure 2E). In summary, these intricate details reveal the 
comprehensive mechanism behind 82:53 activation by the synthetic peptide ВАТ, 


highlighting the specific interactions and critical residues involved in this process. 


BRS3 regulation by MK-5046 

Continuing our exploration of BRS3, we examined its interaction with MK-5046, the 
first BRS3-specific nonpeptide receptor agonist (Figure 3). MK-5046 is characterized 
by distinct functional groups, including a hydroxyl group and a trifluoromethyl group in 
the central segment, an additional trifluoromethyl group accompanied by a 
cyclopropane group and an imidazole group on one side, and a benzene ring and a 
pyrazole group on the other (Figure 3B). The structure of MK-5046-bound BRS3 
reveals that MK-5046 adopts a "V"-shaped conformation, fitting snugly at the bottom 
of the BRS3 pocket. Notably, MK-5046, being smaller in volume (475 À?) than BA1 
(1089 АЗ), occupies a specific sub-pocket that corresponds to the C-terminal HFX(8- 
10)" motif (Figure ЗА). 
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The trifluoromethyl group, connected to the cyclopropane group, projects towards the 
upper portion of TM3, establishing extensive fluorine-based halogen bonding 
interactions with 5124329 212733 and E182* (Figure ЗС). Mutations іп R127??? and 
E182^9? lead to significantly reduced MK-5046-induced signaling (Figures ЗЕ and 58; 
Table S4), emphasizing the critical roles of these residues in MK-5046 binding and 
BRS3 activation. Replacement of К127332 with 0332 in both NMBR and GRPR 
receptors notably decreased the Emax level for mutated BRS3 (Figures 3E and S8; 
Table 54), underlining R127???'s significance in MK-5046's selectivity for BRS3. 


Additionally, the imidazole, benzene, and pyrazole rings of MK-5046 engage in 
extensive п-п interactions with М/284848 and Y291°°° of BRS3 (Figure ЗС). Тһе 
interaction between MK-5046's imidazole group and 2918255 is further reinforced by а 
hydrogen bond with S205*€'? (Figure ЗС). Mutating S205*^? and У291655 results in а 
considerable decrease in MK-5046-induced BRS3 signaling (Figures 3E and S8; Table 
S4), confirming their essential role in this interaction. The central hydroxyl group of 
MK-5046 forms a hydrogen bond with С221542 while the central trifluoromethyl group 
creates robust packing interactions with Н288652 complemented by п-п packing 
between W2849^? and H2885*? (Figures ЗС and 3D). The opposite side of MK-5046, 
featuring a benzene ring and pyrazole group, is accommodated within a hydrophobic 
pocket formed by L96253, М131336, 225546 VV284648, 4319742 and F320’* (Figure 
3D). 


Previous competitive assays involving MK-5046, Peptide-1 (another pan-BnR agonist), 
and Bantag-1 (a BRS3-specific peptide antagonist) suggested that MK-5046 functions 
as an allosteric agonist for BRS3, possibly binding to a distinct pocket from Peptide- 
131, Our structures show that MK-5046's binding site and pose closely resemble the C- 
terminal НЕХ(8-10)841 motif of ВАТ (Figure ЗВ). The upper trifluoromethyl group of MK- 


5046 corresponds to the binding site for BA1's H88, while its central trifluoromethyl 
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group and benzene ring align with the lower hydrophobic cavity, mirroring the FX(9- 
10)? side chains of ВАТ. The pyrazole group of MK-5046 assumes a position akin to 
the carboxyl group of BA1's Х108. In summary, MK-5046's binding pattern shares 
similarities with the C-terminal HFX(8-10)®*' tripeptide of ВАТ, suggesting that MK- 
5046 is an orthosteric agonist. Interestingly, alanine substitution of КЗ16739 enhances 
BRS3 activity induced by MK-5046, showing higher efficacy in our IP1 assay (Figure 
3F and S8; Table S4). This finding further underscores the intricate molecular 


interactions that govern BRS3 regulation and activation by MK-5046. 


Activation of BRS3 by agonists of different types 

Understanding how receptor occupation correlates with activation, particularly in the 
context of peptide and nonpeptide agonists for GPCRs, is crucial. Interestingly, BRS3 
exhibits a unique behavior in this regard. The binding affinities of the peptide agonist 
ВА1 and the nonpeptide agonist МК-5046 are 2.0 nM and 17 пМ1932 respectively, 
indicating that BA1 has a nearly 10-fold higher affinity for the receptor than MK-5046. 
Correspondingly, analysis using dr_sasa reveals that the interface areas for BA1 and 
MK-5046 are 1094.26 А? and 429.02 A?, respectively??. This suggests a direct 
relationship between a larger buried surface area in the structure and a higher binding 


affinity. 


However, the activation potency of these agonists presents an intriguing contrast. In 
our IP1 assays, MK-5046 activates BRS3 with an ECso value of 0.5 nM, demonstrating 
26-fold greater potency in activating Са signaling than ВАТ, which has an ECso of 13.0 
nM. This is notable as most nonpeptide agonists typically exhibit lower potency than 
the natural peptide agonists of their respective receptors in activating downstream 
signaling. Yet, MK-5046 surpasses the peptide ligand BA1 in terms of G, signaling 


activation. 


To further investigate the potency differences between MK-5046-BRS3 and BA1-BRS3, 
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we conducted molecular dynamics (MD) simulations, each spanning 500 ns over three 
replicates. These simulations focused on the ligand-binding domain composed of TM3, 
ТМ5, and TM6, a region known to influence ligand potency**. The simulations revealed 
distinct hydrophobic interaction networks in the extracellular TM3-TM5-TM6 region for 
the two agonists. Specifically, in the MK-5046-BRS3 system, a hydrophobic network 
surrounding the trifluoromethyl group of MK-5046 was present in 62.996 of the 
trajectories, compared to only 23.7% in the BA1 system around F9?^', Representative 
structures (Figure 4A) showed that the trifluoromethyl group in MK-5046 engages in 
electronic interactions with Н2886:52 and facilitates hydrophobic interactions between 
F225546 and М131336. In contrast, the larger Е98А moiety in the ВАТ system interacts 
primarily with М131338, disrupting the hydrophobic contact between 225546 and 
М131336 (Figure 4B). 


Moreover, in the MK-5046 system, the trifluoromethyl group acts as a hydrophobic 
bridge between Н2886:52 and Е225546 stabilizing the activated extracellular TM3-TM5- 
TM6 conformation. This is evidenced by the fact that direct contact between H2885*? 
and Е225° is less frequent (19.9% versus 78.1% in the BA1 system). Consistent with 
this, the H288A mutation significantly diminishes MK-5046's potency (from 0.53 nM to 
126.9 nM), while having a less pronounced effect on BA1's binding (from 11.8 nM to 
68.7 nM) (Figures 2G, 3F, and 4C). In summary, these findings illustrate the complex 
and distinct mechanisms by which different types of agonists, both peptide and 
nonpeptide, activate BRS3. They also highlight the nuanced molecular interactions 
that govern the receptors response to these agonists, contributing to our 


understanding of GPCR activation and signaling. 


Basis for ligands selectivity of bombesin receptors 
In the BnR family, BA1, a synthetic analog, was found to have a high affinity for BRS3, 
as well as for NMBR and GRPR receptors 193558 (Table 51). BRS3 shares 


approximately 50% amino acid sequence identity with NMBR and 4796 with 
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СКРК! 3738. A closer examination of ВАТ binding in BRS3 and GRPR indicates that 
this synthetic peptide uniquely coordinates with all three bombesin receptors. 
Structural alignment of the BA1-GRPR-G, complex (PDB: 7W40)? with our BA1-BRS3- 
G, structure showed a high degree of similarity, with an overall Ca RMSD of 0.75 À 
(Figure 5А). Both the PFQW(2-4)? and HFX(8-10)9^! motifs of BA1 align well between 
BRS3 and GRPR, with notable differences in the side chains of PF2P^!, WA?^', and 
H8®4' (Figure 5B). The key residues in BRS3's upper pocket, including Y 108295, 
Е200ЕС2 H3087?', and Е309732 involved in п-п interactions with PF2P^! and W484", 
are highly conserved in GRPR without substitutions (Figure 5C), suggesting a 
preserved interaction pattern despite conformational differences in these aromatic 
amino acids. However, in BRS3, the H8®4' side chain deviates from TM3, causing a 
shift of the middle AVPA(5-7)?^! motif towards TM6, compared to its position in СЕРЕ 
(Figure 5D). This shift is stabilized by two hydrogen bonds between H8?^! and 5124329 
апа 5205452 in BRS3. In contrast, ВАТ H8 is sandwiched by P??? and P^9*? in СЕРЕ, 
resembling the coordination seen in NMB-NMBR and GRP(14-27)-GRPR structures 
(Figure S7F). The significant sequence homology and specific flexibility, especially in 


the ligand-binding pockets, explain BA1's ability to bind all three bombesin receptors. 


The selectivity of natural bombesin-related peptide agonists for individual bombesin- 
related receptors is an interesting aspect. GRP, the natural ligand for GRPR, has an 
affinity for GRPR over 650-fold higher than for NMBR and more than 15,000-fold higher 
than for BRS3. NMB, the natural ligand for NMBR, shows over 650-fold higher affinity 
for NMBR compared to GRPR and more than 20,000-fold higher than for BRS3 (Table 
S1). Prior studies have clarified the determinants of selectivity between NMB30 and 
GRP(14-27) for ММВК and СКРК 3940, |n our sequence and structure analyses of 
these receptors bound with different peptides, we observed five amino acids in the 
peptide agonist binding pockets that are conserved in both GRPR and NMBR but differ 
in BRS3: T106?65, H107?9^, к127332, 5205:С2, and Н2945 in BRS3 are substituted 
ру 5263, R264, Q332, PEC? and R59, respectively (Figure 51), in NMBR and GRPR. A 
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structure comparison of BA1-BRS3 with NMB-NMBR shows a clash between T6%"® of 
ММВ and 5202522 of BRS3, which is absent in the equivalent interaction in the BRS3 
structure due to a shift of BA1's middle motif towards TM6 (Figure 5E). In the BRS3 
structure, the equivalent V6®"' residue of ВАТ shift away from the 5202542 due to the 
shift of the middle AVPA(5-7)9^! motif towards the TM6. This shift of BA1 middle motif 
may be a structural characteristic of BRS3 specific ligand, which is favored by the 
residue combination of the К127332 and Н294658 substitutions іп BRS3 (Figure 5E). 
Through structural analysis, it is evident that BA1 presents a binding pose of the middle 
motif nearly identical to the natural agonist GRP in the GRPR structures??, located 
close to TM6 (Figure 5F). Interestingly, the NMB middle motif exhibits a closer 
alignment with TM2 (Figure 5F). Consequently, these ligand binding poses of NMB and 
GRP are favored by the conserved residue combination of the Q??? and R$** 
substitutions in NMBR and СЕРКЕ receptors. Thus, the Q??? апа К substitutions 
contribute to the structural characteristics of ligand specific for individual Bn-related 
receptors, as suggested in the prior mutation studies^'. For ММВ, it has been reported 
that T6"V? forms a hydrogen bond with the main carboxyl group of L3“™®, playing a 
pivotal role in ММВ binding and selectivity by NMBR*?, The clash of the T6^V? here 


possibly contributes to the weak binding of the NMB by BRS3. 


In addition, structure comparison of the BA1-BRS3 with our previously reported 
GRP(14-27)-GRPR structure (PDB code: 8H0Q) reveals that the T106?9? substitution 
in BRS3 induces an inward push of the side chain of W113E°™ (Figure 5G). This 
movement is further stabilized by the packing between the side chain of W1135°" and 
the H1072 substitution in BRS3, along with the conserved disulfide bond formed by 
С120325 and C203°°!? (Figure 5G). Together with the main chain carboxyl group of 
C201-%X2 the inward shift of W113ECH subsequently results іп an inward movement of 
the EGW(111-113)***! of BRS3, causing a clash with the crucial residue H39"P of 
GPR(14-27) (Figure 5H). When we compare the BA1-BRS3 with the other GRP-GRPR 
structure (PDB code: 7W3Z), a clash between the H39"* and the carboxyl group of 
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E2015? can also be observed (Figure 576), despite the H3°"? here adopt a rotemer 
conformation due to the similar ECL1 positions between the BRS3 and GRPR 
structures. Interestingly, the H3°"? was also reported to accommodate the featured 
ECL2 of GRPR, and contribute to the weak binding of GRP to NMBR???, suggesting 
its pivotal role of selectivity for GPR(14-27) in BnR family. Overall, these structural 
insights highlight the unique features of BRS3 that contribute to its low affinity for 


bombesin-related natural ligands like NMB and GRP. 


Common activation features of bombesin receptors subfamily 

In examining the significant structural similarities between the BA1- and MK-5046- 
bound 8:53 complexes, we chose the BA1-BRS3 complex for a detailed analysis of 
BRS3 activation. Comparing the active BA1-bound BRS3-G4 complex with both the 
active structure of GRPR bound to BA1 (PDB: 7VV40) and the inactive structure of 
GRPR bound to a synthetic antagonist, PD176252 (PDB: 7W41), we observe that 
BRS3 exhibits classic activation features typical of Class A GPCRs in response to ВАТ. 
Notably, the cytoplasmic end of TM6 in 8:53 displays a pronounced outward 
movement while the cytoplasmic part of TM7 shifts inward. This rearrangement 
accommodates the C-terminal a5 helix of the Ga, subunit, a key feature of Class A 
GPCR activation (Figures 6A-6C). Additionally, the extracellular end of TM1 in the BA1- 
bound BRS3 shifts 5.3 A towards TM7, similar to movements observed in NMBR and 
GRPR structures (Figures 6B and S7H). 


In comparison to PD176252, BA1 binds deeper within the orthosteric pocket, 
surrounded by ТМЗ, TM6, апа ТМ7, and forms hydrophobic contacts with W284°*8 
(Figure 6D). BA1's binding causes a significant deflection of the VV284648 indole ring, 
initiating a cascade of conformational changes linked to receptor activation’? (Figure 
6D). This process mirrors the classic activation mechanism seen in Class A GPCRs, 
where the displacement of М/284848 exerts a downward force оп Е28084 in the 


conserved Р950у/340р644 core triad, leading to changes in the D9^9?R350Y$5! motif 
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(Figures 6D-6F). Additionally, the conserved N’49P7°°xxY753 motif is replaced with a 
N’49P7-5°xxL7-53 motif іп BRS3. In the structure, BA1's X10*^! carboxyl group forms a 
hydrogen bond with R3167??, contributing to activating the М749Р759хх 753 motif in 
BRS3. Typically, the conserved Y^? in this motif undergoes a significant rotation 
towards TM3, leading to TM7’s inward movement*?. However, іп BRS3, the leucine 
substitution at Y’°? undergoes a lateral shift towards TM1 and TM2, highlighting a 


unique yet conserved activation mechanism shared among bombesin receptors. 


During BRS3 activation, the P°5°V34°F&44 and N749P750xx| 755 motifs, which link the 
ligand-binding pocket to the G protein-coupling interface, undergo rearrangements. 
These involve movements of F844 and L7:53, triggering significant outward shifts of TM5 
and TM6 (3.4 A at К°; 8.1 A at Ке), thereby facilitating the engagement of the Сод 
subunit's C-terminus. This feature is consistent with activation mechanisms observed 
across GPCRs (Figure 6C; Table S5). These structural insights from our BRS3 studies 
enhance our understanding of the activation mechanisms within the bombesin receptor 


subfamily, demonstrating how specific molecular changes lead to receptor activation. 


Comparison with the AlphaFold2-predicted structure 

We compared our cryo-EM resolved structures with the structure predicted by 
AlphaFold2, specifically focusing on the receptor domain (Figure 571). Given the 
notable structural similarities between the BA1- and MK-5046-bound BRS3 complexes, 
we selected the BA1-BRS3 complex for further analysis. The resolved BA1-BRS3 
structure and the AlphaFold2-predicted structure exhibited somewhat dissimilar overall 
arrangements, with a Ca RMSD of 1.5 A. Despite AlphaFold2 capturing the general 
backbone contours of the receptor, substantial discrepancies existed between the 
predicted models and experimental structures regarding the assembly of the 
extracellular and transmembrane domains. Compared to the BA1-BRS3 structure, the 
AlphaFold2 structure represents an inactive conformation, with the intracellular part of 


TM6 tilting outward and the intracellular part of TM6 shifting inward. Specifically, the 
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key residues іп the ОВР in the AlphaFold2 8:53 structure, including Н107264 5124329, 
R127?*?, W284°*8 (Figure S7J), adopt distinct positions that cannot further favor the 
BA1 binding. In conclusion, while the AlphaFold2-predicted model serves as a valuable 
reference point, obtaining additional empirical structural information is imperative to 


accurately guide drug design. 


Discussion 


Our study on BRS3, an orphan member of the bombesin receptor subfamily, offers 
significant insights into the complex dynamics and mechanisms of GPCRs. As orphan 
receptors, entities like BRS3, without identified endogenous ligands, present 
challenges in understanding their roles, particularly in physiological and pathological 


contexts like metabolic diseases and cancer^^^$6. 


The structural analysis of BRS3 in its apo form and when bound to the synthetic analog 
BA1 and the nonpeptide agonist MK-5046 reveals a high degree of structural 
conservation (Figure S7A). This similarity, despite the evolutionary loss of binding to 
natural ligands like NMB and GRP in placental mammals, underscores the evolutionary 
adaptability of BRS3*’. Unlike other orphan receptors that drive constitutive activation 
through their ECL2s, such as GPR12 and GPR2148°°, BRS3 does not exhibit this 


feature, indicating unique structural and functional characteristics. 


Our investigation reveals unique aspects of ligand recognition and activation in BRS3. 
The complex interplay of residues within the receptor's binding site, including those 
contributing to agonist-induced activation, highlights the receptor's intricate regulation 
mechanisms. These findings provide a foundation for the future exploration of BRS3 
as a potential therapeutic target, particularly in diseases linked to dysfunctional GPCR 


signaling. 


In the bombesin receptor family, the synthetic analog BA1 exhibits high affinity for 


BRS3, NMBR and GRPR. Our structural analysis comparing the BA1 complexes with 
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8:53 and СЕРЕ shows that the overall binding mode of BA1 to the two receptors is 
highly similar despite some differences in their amino acid side chains. This is largely 
attributable to the high homology of key residues in the ligand binding pockets of these 
receptors. These findings explain BA1's ability to bind all members of the BnR family. 
On the other hand, natural peptides like NMB and GRP demonstrate high selectivity 
for individual BnR receptors. Our sequence and structural comparisons of apo and 
peptide-bound BRS3 identified five key residues that differ from NMBR and СЕРЕ, 
which may account for its low affinity for NMB and GRP. This provides evidence to 
understand the loss of binding to these natural peptides by BRS3. Our study offers 
insights into the structural basis of ligand selectivity, aiding efforts to design BRS3 


peptide ligands with high selectivity. 


In addition to these findings, our research aligns with the classical activation features 
observed in Class A GPCRs, including the outward movement of TM6 and the inward 
shift of TM7 in response to BA1 binding. This conformational change, particularly the 
deflection of the W284°*® indole ring, echoes activation mechanisms seen in other 
GPCRs. Such structural insights into BRS3, in tandem with its evolutionary trajectory, 
contribute to a deeper understanding of its functions and potential therapeutic 


implications. 


Our comprehensive study, therefore, not only unravels the unique structural and 
functional attributes of BRS3 but also enhances the overall understanding of GPCR 
activation mechanisms. This knowledge is crucial for developing targeted interventions 
for diseases associated with GPCR dysfunction. Future research is needed to delve 
deeper into BRS3's self-activation mechanism and explore its potential as a 


therapeutic target within the bombesin receptor family. 
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675 
676 Figure 1. Cryo-EM structures of BRS3-G4 complexes, both in the absence (apo- 


677 form) and presence of the ligands ВАТ and MK-5046. 

678 (А) Schematic illustration of the BRS3 function in the regulation of metabolism 

679 pathways. Briefly, activation of BRS3 neurons in the hypothalamus increased 

680 sympathetic nervous tone, elevating energy expenditure, body temperature, heart 
681 rate, and blood pressure via brown adipose tissue, heart, and arteries, respectively. 


682 Апа agonists activate BRS3 on pancreatic В cells, leading to an increase in glucose- 
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stimulated insulin secretion. 

(B and C) Orthogonal views of the density map (B) and model (C) for the apo BRS3- 
Сч complex. The receptor BRS3 is displayed in forest green. The heterotrimeric бұ 
proteins are colored by subunits: Gag, yellow; GB, cornflower blue; Gy, deep salmon; 
scFv 16, gray. 

(D and E) Orthogonal views of the density map (D) and model (E) of the BA1-BRS3- 
G, complex. BA1 is shown in light salmon; the receptor BRS3 is displayed in 
aquamarine. The heterotrimeric ба proteins are colored by subunits: Gag, yellow; GB, 
cornflower blue; Gy, deep salmon. 

(F and G) Orthogonal views of the density map (F) and model (G) of the MK-5046- 
BRS3-G, complex. MK-5046 is shown in orange; the receptor BRS3 is displayed in 
violet. The heterotrimeric Са proteins are colored by subunits: Саг, yellow; СВ, 


cornflower blue; Gy, deep salmon; scFv16, gray. 
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Figure 2. The BA1-binding pocket in BRS3. 

(A) The sequence of BA1 peptide. The amino acid numbering of BA1 is determined 
relative to NMB. See also Table S3. 

(B) Side views of the ligands-binding pocket in bombesin receptors family. 

(C-E) Detailed interactions of BA1 with BRS3. Hydrogen bonds are depicted as black 
dashed lines. 

(F and G) Effects of mutations in the BA1-binding pocket in IP1 assay. Data are 
presented as mean + S.E.M. of at least three independent experiments. Source 


data are available online. 
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740 Figure 3. The MK-5046-binding pocket in BRS3. 


741 (А) The chemical structure of MK-5046 and detailed interactions between МК-5046 


742 and BRS3. 
26 


743 


744 


745 


746 


747 


748 


749 


750 


751 


752 
753 
754 
755 
756 
757 
758 
759 
760 
761 
762 
763 
764 
765 
766 
767 
768 
769 
770 
771 
772 
773 
774 
775 
776 
777 


(В) The structural comparison of BA1-BRS3 and MK-5046-BRS3 complexes. Тһе 
binding site and pose of MK-5046 closely resemble the C-terminal dumbbell end 
HFX(8-10)®4' motif of ВАТ. 

(C and D) Detailed interactions of MK-5046 with BRS3. Hydrogen bonds are depicted 
as black dashed lines. 

(E and F) Effects of mutations in the MK-5046-binding pocket in IP1 assay. Data are 
presented as mean + S.E.M. of at least three independent experiments. Source data 


are available online. 
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Figure 4. Activation of BRS3 by agonists of different types 

(A and B) The representative structures revealed the difference in the MK-5046- 
BRS3 (A) and BA1-BRS3 (B) complexes. 

(C) Effects of H288°*°2A on the potency in IP-1 assay. Data were shown as ApECso + 
S.E.M. from three independent experiments, which performed in triplicates, with total 
repeats of nine for each data point. Statistical differences between wild type (WT) 
and mutants were determined by two-sided, one-way analysis of variance (ANOVA) 


with Tukey test. *, P < 0.05; **, P < 0.01; ***, ჩ < 0.001; ns, no significant difference. 
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Figure 5. Molecular basis of ligands selectivity of bombesin receptors. 

(A) A structural comparison of BA1-BRS3 and BA1-GRPR complexes. The ВАТ in 
BA1-BRS3 complex is shown in light salmon; the receptor BRS3 is displayed in 
aquamarine; the ВАТ in BA1-GRPR complex is shown in dodger blue; the receptor 
GRPR is displayed in burly wood. 

(B) A structural comparison of BA1 conformation in BA1-BRS3 and BA1-GRPR 
complexes. 

(С and D) Detailed interactions between PF2P^', WA4P^'(C), and Н ВМА (D) with 


residues in СЕРЕ and BRS3. The shift of middle AVPA(5-7)9^' motif is represented 
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by an arrow. 

(E) Detailed interactions between Т6 98, V6ERP, and V6B^ with residues іп NMBR, 
СЕРЕ and BRS3. The shift of Т6УМВ is represented by an arrow. 

(Е) Detailed interactions between | ЗУУВ, H3SFP, and ОЗ8А with residues in NMBR, 
GRPR and BRS3. The middle motif of BA1 and GRP exhibits a closer alignment with 
TM6, and the counterpart of NMB closed to TM2. 

(С) The inward push of the side chain of W113ECH by Т106263 substitution and the 
packing between the side chain of №1138 and the H1072 substitution in BRS3. 
(Н) The clash between H3°"? of GPR(14-27) апа EGW(111-113)*€*! of BRS3 is 


represented by a circle. 
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Figure 6. Molecular basis for activation of BRS3. 

(A-C) Structural alignment of two active BRS3 in the presence of the ligands BA1 and 
MK-5046, active BA1-GRPR (PDB code: 7W40) and inactive GRPR (PDB code: 
7\41). Side view (A); extracellular view (B); cytoplasmic view (C). The movement 
directions of TM6 апа TM7 in BRS3 relative to inactive GRPR are highlighted as red 
arrows. Inactive GRPR, active GRPR, active BRS3 in the presence of the ligands 
BA1 and MK-5046 are colored in gray, burly wood, aquamarine and violet, 
respectively. 

(D-F) Conformational changes of the conserved 'micro-switches' upon receptor 
activation, including toggle switch (D), PV(I)F (E), and DRY(F) motifs. The 
conformational changes of residue side chains are shown as red arrows upon 


receptor activation. The complex structures were aligned by the receptors. 
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Methods 
Constructs 


Human BRS3(residues 1-361) was cloned into pFastBac with an N-terminal 
haemagglutinin (HA) signal peptide followed by thermostabilized apocytochrome 
bse2RIL (BRIL)'. LgBiT was also inserted at the C terminus of the BRS3 using 
homologous recombination (CloneExpress One Step Cloning Kit, Vazyme). After 
LgBiT, a tobacco etch virus (TEV) protease cleavage site and a double maltose 
binding protein (MBP) tag were added to 8553 constructs for improved expression 
and purification. The engineered Gag construct was generated based on mini-Gsq7 12, 
which carries two dominant-negative mutations (G203A and A326S)? to decrease 
nucleotide binding affinity and facilitate GaBy complex stability. The N-terminal 1-18 
amino acids and the a-helical domain of the mini-Gq71 were substituted with the 
corresponding sequences of the human Gai for binding to the antibody fragments 
scFv16*°. Rat GB1 was connected with a C-terminal HiBiT by а 15 residues linker. 
The engineered Ga, ОВ1-НІВІТ and bovine Gy2 were cloned independently into 
pFastBac vectors (Invitrogen). The antibody fragment scFv16 was cloned into a 


modified pFastBac vector containing an N-terminal GP67 secretion signal peptide. 


Expression and purification of Nb35 

Nanobody-35 (Nb35) with a N-terminal pelB signal peptide and a C-terminal Hise tag 
was expressed in the periplasm of Escherichia coli BL21(DE3) bacteria (МЕВ). 
Cultures were grown at 37 °C in Luria-Bertani media containing 50 ug/mL ampicillin 
to an ODeoo of 1.0 and induced with 0.1mM IPTG at 28°C, 180 r.p.m. for another 8 ი. 
Cells were harvested by centrifugation (5316 x g, 30 min) and lysed in ice-cold buffer 
(20mM HEPES pH 7.4, 500mM NaCl), then centrifuged to remove cell debris. No35 
was first purified by nickel affinity chromatography, followed by size-exclusion 
chromatography using a HiLoad 16/600 Superdex 200 prep grade column with a 
buffer of Zo0mM HEPES pH 7.4 and 100mM NaCl. Peak fractions were concentrated 
to 2mg/mL with 15% glycerol and kept frozen at -80 °C for later use. Purified protein 


quality was assessed by SDS-PAGE. 
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Expression and purification of complexes 


High Five cells (Expression systems) were cultured in ESF921 serum-free medium 
(Expression Systems) and co-infected with baculoviruses expressing the receptor, 
Саа, СВ1, Gy2, and scFv16' at a 1:1:1:1:1 ratio for 48 h at 27 °С using Bac-to-Bac 
baculovirus system. The cell pellets were lysed by dounce homogenization in a buffer 
containing 20 mM HEPES pH 7.4, 100 mM NaCl, 10 mM MgCle, 0.1 mM TCEP 
(Sigma-Aldrich), 1096 glycerol, and EDTA-free protease inhibitor cocktail (TargetMol). 
The supernatant was isolated by centrifugation at 65,000 x g for 40 minutes to collect 
the membranes. The washed membranes were resuspended in buffer containing 20 
mM HEPES pH 7.4, 100 mM NaCl, 10 mM MgCle, 10% glycerol, agonists (10 uM 
BA1, and 10 uM MK-5046, Synpeptide and MedChemExpress, respectively), 25 
mU/mL apyrase (Sigma-Aldrich), 20ug/mL Nb35, 0.1 mM TCEP and EDTA-free 
protease inhibitor cocktail, and incubated at 4°C for 8 hours. After incubation, cell 
membranes were solubilized with 1% (w/v) n-dodecyl-8-d-maltopyranoside (DDM, 
Anatrace) and 0.296 (w/v) cholesteryl hemisuccinate TRIS salt (CHS, Anatrace) at 
4°C for 2 hours. The supernatant was collected by centrifugation at 65,000 x g for 35 
minutes and incubated with dextrin resin (Dextrin Beads 6FF, Smart Life Sciences) at 
4'C for 5 hours. The resin was washed, eluted, and the protein complexes were 
further purified by size-exclusion chromatography on a Superdex 200 increase 
10/300 GL column (GE Healthcare) in buffer containing 20 mM HEPES pH 7.4, 100 
mM NaCl, 2 mM MgCl», 0.1 mM ТСЕР, agonists (10 uM BA1, 10 uM MK5046, and 1 
mM oridonin), 0.00075%(w/v) LMNG, 0.00025% (w/v) GDN and 0.0002% (w/v) CHS. 
Fractions containing monomeric protein complex were collected and evaluated by 
SDS-PAGE (Figure S2), and then concentrated by 30-fold for cryo-electron 


microscopy experiments. 


Cryo-EM data collection 
Cryo-EM grids were prepared with the Vitrobot Mark IV plunger (ЕЕ!) set to 8 "С and 


100% humidity. Three-microliter of the BRS3-G, protein complex in apo state was 
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applied to the glow discharged copper R1.2/1.3 holey carbon grids. The sample was 
incubated for 105 on the grids before blotting for 3.5 5 (double-sided, blot force 1) and 
flash-frozen in liquid ethane immediately. The same condition was used for the sample 


BA1-BRS3-G, complex, and MK5046-BRS3-G, complex. 


For BRS3-G, complex dataset, 3,098 movies were collected on a Titan Krios equipped 
with a Gatan K3 direct electron detection device at 300 kV with a magnification of 
64,000, corresponding to a pixel size 1.08 A. Image acquisition was performed with 
EPU Software (FEI Eindhoven, Netherlands). We collected a total of 36 frames 


accumulating to a total dose of 50 е: А? over 2.5 s exposure. 


For BA1-BRS3-G, complex dataset, 5,016 movies were collected on a Titan Krios 
equipped with a Gatan КЗ direct electron detection device at 300 kV with a 
magnification of 105,000, corresponding to a pixel size 0.824 A. Image acauisition was 
performed with EPU Software (FEI Eindhoven, Netherlands). We collected a total of 


36 frames accumulating to a total dose of 50 е: А? over 2.5 s exposure. 


For MK5046-BRS3-G, complex dataset, 8,738 movies were collected on a Titan Krios 
equipped with a Falcon4 direct electron detection device at 300 kV with a magnification 
of 165,000, corresponding to a pixel size 0.73 A. Image acquisition was performed with 
ЕРІ) Software (FEI Eindhoven, Netherlands). We collected a total dose of 50 е: А? 
over 2.5 s exposure on each EER format movie®. Each movie was divided into 36 


frames during motion correction. 


Cryo-EM image processing 
MotionCor2 was used to perform the frame-based motion-correction algorithm to 
generate drift-corrected micrograph for further processing, and CTFFIND4 provided 


the estimation of the contrast transfer function (CTF) parameters?'°. 
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For BRS3-G, complex dataset, approximately 2,000 particles were manually picked 
and two-dimensional (2D) classes were calculated and used as references for 
automatic picking. All subsequent steps including particle picking and extraction, 2D 
classification, three-dimensional (3D) classification, 3D refinement, CTF refinement, 
Bayesian polishing, post processing and local resolution estimation were performed 


using Relion3.0!!. 


A total of 3,294,877 particles were extracted from the cryo-EM micrographs and 
followed by two rounds of reference-free 2D classification, yielding 3,270,779 particles 
after clearance. Three rounds of 3D classification were used to separate out 233,381 
particles, which were refined to a structure at 3.93 À global resolution. After CTF 
refinement, polishing, and postprocessing, the particles were reconstituted to a 3.29 À 


structure. (Figure S3; Table S2). 


For BA1-BRS3-G, complex dataset, the subsequent steps were performed with 
Сгуо5РАНС?2. Blob-pick was used for the particle picking. 6,245,406 particles were 
automatic picked and extracted from 5,016 cryo-EM micrographs. After 2 rounds of 2D 
classification, 799,895 particles were selected from 2,116,883 particles and subjected 
to generate 5 "good" references and 1,040,145 particles were selected and subjected 
to generate 5 "bad" references. Two rounds of heterogeneous refinement using "good" 
and “Бад” references yielded 675,060 particles after clearance. We then continued the 
processing in CryoSPARC. After homogeneous refinement and non-uniform 


refinement, the particles were reconstituted to a 2.90 À structure (Figure S4; Table S2). 


For datasets of MK-5046-BRS3-G, complex, the similar strategy was used. After non- 
uniform refinement, 196,746 particles were reconstituted to 2.93 À structure, which 
corresponding to MK-5046-BRS3-G, complex (Figure S5; Table S2). 


Model building 
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BRSS structure predicted from Alphafold2 was used as the starting reference models 
for receptors building'?. Structures of Сол, СВ, Gy, the NB35 nanobody, and scFv16 
derived from PDB entry 8HO0Q'^ were rigid body fit into the density. All models were 
fitted into the EM density map using UCSF Chimera’® followed by iterative rounds of 
manual adjustment and automated rebuilding in СООТ' and РНЕМІХ"”, respectively. 
The models were adjusted in ISOLDE'? followed by refinement in PHENIX. The final 
model statistics were validated using Comprehensive validation (cryo-EM) in 
РНЕМІХ"” and provided іп the Table S2. All structural figures were prepared using 
Сһітега!5, Chimera X'?, and PyMOL (Schrödinger, LLC.). 


Function assay 

AD293 cells (Agilent) were cultured in high glucose DMEM medium (GE healthcare) 
supplemented with 1096 (v/v) fetal bovine serum (FBS, Gemini) and 196 
penicillin/streptomycin at 37°C in 5% CO» incubator. Inositol phosphate 1 (IP1) 
production was measured using the IP-One HTRF kit (Cisbio, 621РАРЕЈ)2°. Briefly, 
cells were seeded onto 12-well plates 16 hours before transfection. Cells were then 
transiently with different BRS3 constructs using FuGENE HD transfection reagent. 
After 24 hours, cells were harvested and resuspended at 7 x 10? cells/mL in IP1 
stimulation buffer. Cells were then plated onto 384-well assay plates at 4,900 cells/7 
uL/well. 7 uL IP1 Stimulation Buffer 2 containing ligand was added and incubated for 
1 hour at 37 °C. Intracellular ІРІ was measured using the IP-One HTRF kit and an 
EnVision multiplate reader according to the manufacturer's protocol. The HTRF ratio 
was converted to a response (%) using: response (Фо) = ratio of sample/WT x100. 
Data presented are mean + S.E.M. of at least three independently biological 


experiments. 


Cell-surface expression assay 
Cell-surface expression for each 8553 mutant was monitored by a fluorescence- 


activated cell sorting (FACS). The mutants were cloned into pcDNA6.0 vector 
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(Invitrogen) with a N-terminal FLAG tag. Cell were seeded and transfected as 
described for the functional assay. After 24 hours of transfection, cells were washed 
with PBS and detached with 0.2% (w/v) EDTA in PBS. The expressed cells were 
then incubated with monoclonal anti-FLAG M2-FITC (Sigma-Aldrich) at 1:100 dilution 
for 15 minutes at 4 °C, followed by a 9-fold excess of PBS. Cells were resuspended 
and fluorescence intensity was quantified using a BD Accuri C6 flow cytometer (BD 
Biosciences) with 488 nm excitation and 519 nm emission. FACS data were 


analyzed with BD Accuri C6 software 1.0.264.21 and normalized to wild-type 886553. 


Computational methods 

The simulation systems come from MK-5046-BRS3-G protein complex and ВА1- 
BRS3-G protein complex. G proteins were removed before simulations. Protonation 
states of residues were determined using Propka3 software?'. The CHARMM-GUI 
platform facilitated embedding of these structures into a 75x75 A POPC lipid bilayer”. 
The lipid bilayer was then surrounded by a 15 À aqueous layer. Systems were then 
adjusted to a 0.15 mol/L NaCl concentration, supplemented with counterions. We 
applied the CHARMM36m force field for amino acids and lipids, while CHARMM 
general force field is applied for ligands?*?^. The systems underwent a 7-step 
equilibration process, with gradual minimization and relaxation of constraints, as 
outlined by CHARMM-GUI. We conducted 3x500 ns independent production runs for 
each system using pmemd.cuda in Amber20? under the NPT ensemble at 303.15 K 
and 1 atm. Long-range electrostatics were managed via the Particle Mesh Ewald 
method, while short-range electrostatic and van der Waals interactions used a 12 À 
cutoff, with a smooth transition between 10 and 12 À. The hydrophobic network was 
evaluated using “nativecontact” command in CPPTRAJ. The interacting interface area 


was calculated using dr. sasa?$. 
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1128 
1129 Figure S1. Sequence alignment of BRS3, NMBR, and GRPR. Secondary structure 


1130 elements are annotated underneath the sequences based on the structure of the 
1131  BRS3-G4 complex. The ligand-binding pockets of NMB-NMBR, GRP(14-27)-GRPR, 
1132 BA1-BRS3, and MK-5046-BRS3 is displayed in different colors, purple, orange, 
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1133 green, and violet, respectively. Additionally, the non-conserved residues including 


1134 5263 R264 (0332 р4552 апа В, is circled in dotted box. 
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Figure S2. Purification and characterization of BA1-BRS3-G, complex and MK- 
5046-BRS3-G, complex. 

(A) Schematic diagram of the protein engineering of BRS3, engineered Gag, GB, Gy, 
scFv16, and Nb35 used in this study. 

(B) Gel filtration (Superdex 200 Increase 10/300 column, GE Healthcare) profile of the 
apo BRS3-G, complex. The black arrow indicates the sharp peak for аро BRS3-Gq 
complex. 

(C) Coomassie-stained SDS-PAGE analysis of the purified apo BRS3-G4 complex, 
showing balanced ratios for each subunit. 

(D) Gel filtration (Superdex 200 Increase 10/300 column, GE Healthcare) profile of the 
BA1-BRS3-G4 complex. The black arrow indicates the sharp peak for BA1-BRS3—G, 
complex. 

(E) Coomassie-stained SDS-PAGE analysis of the purified BA1-BRS3-G4 complex, 
showing balanced ratios for each subunit. 

(F) Gel filtration (Superdex 200 Increase 10/300 column, GE Healthcare) profile of the 
MK-5046-BRS3-G4 complex. The black arrow indicates the sharp peak for MK-5046- 
BRS3-G, complex. 

(G) Coomassie-stained SDS-PAGE analysis of the purified MK-5046-BRS3-G, 


complex, showing balanced ratios for certain subunits, except for the Gy subunit. 
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Figure S3. Cryo-EM data processing of apo BRS3-G, complex. 


(A) Flowchart of computational sorting of cryo-EM data. 


(В) A representative cryo-EM micrograph of apo BRS3—G, complex with 50 nm scale 


bar included as a size reference. 


(С) Twelve representative reference-free 2D cryo-EM class averages. Scale bar, 5 пт. 
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(О) ‘Gold-standard’ Fourier shell correlation curve of the reconstruction. The resolution 
was reported at 3.30 A using the Fourier shell cutoff at 0.143. 


(E) Local resolution map of apo BRS3—G, complex. 
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1199 Figure S4. Cryo-EM data processing of BA1-BRS3-G, complex. 

1200  (A)Flowchart of computational sorting of cryo-EM data. 

1201  (B)Arepresentative cryo-EM micrograph of BA1-BRS3-G4 complex with 50 nm scale 
1202 bar included as a size reference. 

1203 (C) Twelve representative reference-free 2D cryo-EM class averages. Scale bar, 5 nm. 


1204  (D)'Gold-standard' Fourier shell correlation curve of the reconstruction. The resolution 
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was reported at 2.90 A using the Fourier shell cutoff at 0.143. 


(Е) Local resolution map of BA1-BRS3—G, complex. 
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1230 Figure 55. Cryo-EM data processing of MK-5046-BRS3-G, complex. 
1231 (А) Flowchart of computational sorting of cryo-EM data. 
1232  (B)Arepresentative cryo-EM micrograph of MK-5046-BRS3-G, complex with 50 nm 
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scale bar included as a size reference. 

(C) Twelve representative reference-free 2D cryo-EM class averages. Scale bar, 5 nm. 
(D) 'Gold-standard' Fourier shell correlation curve of the reconstruction. The resolution 
as reported at 2.93 À using the Fourier shell cutoff at 0.143. 

(E) Local resolution map of MK-5046-BRS3-G, complex. 
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Figure 56 
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1260 
1261 


1262 Figure 56. Cryo-EM density maps of TM1-7 in receptors, aH5 and aHN in Са 
1263 proteins, and ligands in apo BRS3-G, complex structure(A), BA1-BRS3-G,4 
1264 complex structure (В), and MK-5046-BRS3-G, complex structure (С). 
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Figure 57 
A 


apo BRS3 
5 ТМ6 BA1-BRS3 
MK-5046-BRS3 Extracellular view 


т. / 
H H ae 4 
: ot A 
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назхунвумвунвске 
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BA1-BRS3 
тмв NMB-NMBR 
a GRP(14-27)-GRPR 


BRS3 in BA1-BRS3 


NMBR in NMB-NMBR nd 
GRPR in GRP(14-27)-GRPR BA1-BRS3 
C | : BRS3 in AlphaFold2 


1275 
1276 
1277 Figure S7. Structural comparison of bombesin receptors family. 

1278 (A-C) Different views of the structural superposition of apo-state BRS3 (forest green), 
1279  BA1-bound BRS3 (medium aquamarine), MK-5046-bound BRS3 (dark khaki). 

1280 (О and E) Side view of the structural superposition of BA1-BRS3, NMB-NMBR, and 


1281 GRP(14-27)-GRPR(D), and three ligands(E). 
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(Е) The H8? residue forms hydrogen bonds with 5124329 and 5205882, differing from 
its analogs in NMB and GRP(14-27). 

(С) The structural superposition of BRS3 in BA1-BRS3, NMBR in ММВ-ММВК, GRPR 
in GRP(14-27)-GRPR, and GRPR in PD176252-GRPR. 

(H and I) The structural comparison of BRS3 in BA1-BRS3 and AlphaFold2 website(H) 
and the clashes between ВАТ апа ОВР in the AlphaFold2 BRS3 structure(l). 
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Figure 58 
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1312 
1313 Figure S8. ІР1 response curves of ВАТ and MK-5046. BRS3 mutations on agonists 


1314  BA1 (A-E) and МК-5046 (Е-1) induced signaling effects. Dose-response curves for 
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ligands were performed in transiently transfected AD293 cells. The results are 
expressed as percentage of the maximal effect observed in BRS3 wild type. Wild-type 
and mutated BRS3 data are shown by different colors as figures. Each point represents 


the mean + S.E.M. of 3 to 7 independent experiments with triplicate determinations. 
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Figure 59 
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Figure 59. Sequence alignment of BRS3 of different species including 


placental mammals (human, macaca mulatta, rabbit, rat, and mouse) and non- 


placental vertebrates (chicken, turtle, snake, and frog). 
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Table S1. The affinity(pKi) data of Bn-related ligands for human bombesin receptor 


family. 
pK? 
NMB GRP(14-27) ВАТ MK-5046 
hNMBR 8.1-10.3 6.8-6.9 7.8-9.7 «5.0 
hGRPR 7.7-7.8 9.7-9.8 7.8-8.0 <5.0 
hBRS3 5.7 <5.0 8.1-8.4 7.7-8.4 


“Тһе data are from the website “https://www.guidetopharmacology.org/”. 
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Data collection and 


processing 
Detector 


Magnification 

Voltage (kV) 

Electron exposure (e7/A?) 

Defocus range (ит) 

Pixel size (A) 

Symmetry imposed 

Initial particle projections (по.) 

Final particle projections (no.) 

Map resolution (A) 

Map resolution range (A) 
FSC threshold 


Model Refinement 
Refinement package 


Real or reciprocal space Real space Real space Real space 
Model-Map CC (mask) 0.43 0.63 0.62 
Model resolution (A) 3.30 2.87 3.70 

FSC threshold 0.143 0.143 0.5 
8 factors (A2, min/max/mean value) 

Protein residues 22.81/126.95/65.54 20.00/126.95/64.51 22.81/155.34/71.34 

Ligands -H- 20.00/20.00/20.00 20.00/20.00/20.00 
Model composition 
Non-hydrogen atoms 9,067 7,380 10,071 
Protein residues 1,153 930 1,281 
R.m.s. deviations 

Bond lengths (À) 0.001 0.003 0.003 

Bond angles (°) 0.386 0.649 0.924 
Validation 
MolProbity score 1.24 1.26 1.19 
Clashscore 4.70 5.01 4.09 
Rotamer outliers (96) 0.10 0.25 0.09 
Ramachandran plot 

Favored (96) 98.59 98.47 98.42 

Allowed (96) 1.41 1.53 1.58 

Disallowed (96) 0.00 0.00 0.00 
Data availability 
EMDB entry EMD-38927 EMD-38928 EMD-38929 

PDB entry 85 BSS 


apo BRS3-G,- 


scFv16 complex 


кз 
64,000 
300 
50 
-1.0—-3.0 
1.08 
C1 
3,294,877 
233,381 
3.3 
2.60-4.40 
0.143 


PHENIX-1.17.1-3660 
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BA1-bound 
ВК53-Са 
complex 


-1.0~-3.0 
0.824 
C1 
6,245,406 
675,060 
2.90 
2.40-4.20 
0.143 


PHENIX-1.17.1-3660 


Table S2. Cryo-EM data collection, model refinement and validation statistics. 


MK-5046-bound 
BRS3-Go- 
scFv16-Nb35 
complex 


Falcon 4 
165,000 
300 
50 
-1.0~-3.0 
0.73 
C1 
3,416,175 
196,746 
2.93 
2.40-4.20 
0.143 


PHENIX-1.17.1-3660 
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Table S3. Sequence of Bn-related peptides used or mentioned in this study. 


position relative to Bn 


Peptide name N-Terminus 1 уй 8 8 3 9 40 11 12 13 44 
Endogenous 
Bn(Bombesin) E О R L см QW A V б Н 1 мм, 
GRP27(GRP(1-27)) VPLPAGGGTVLTK M Y P К G N H W A V С Н Г ММН, 
GRP (14-27)a(GRP (14-27)) M Y P R С N H W A V G Н L ММН. 
GRP (14-27)b(GRP (14-27)) M Y P R G N H W A V G H L M 
NMB(Neuromedin B) G М М A T G Н F MNW 
NMB30(Neuromedin B (1-30)) LSWDLPEPRSRASKIR V H R R G N L W A T G H F ММН, 
NMC(Neuromedin C) G N H W A V G H L MNH 
Exgenous 

6 13 
Шея OE @ МА V АНЕ X 
Phe “,Ме "]bombesin-(6-14)) 
1 2 3 4 8 8 T 8 9 10 


position relative to NMB or MMC 
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Table S4. Ligand binding affinities and expression levels of WT and mutated 
BRS3. The wild type (WT) and mutants of BRS3 discussed in this manuscript were 
individually analyzed. The affinities are derived from at least 3 independent 
experiments using IP1 function assay. The expression level of mutant BRS3 were 
normalized to wild-type BRS3 as 100%, respectively. Each data point represents mean 
+standard error of the mean (S.E.M.). All data were analyzed by two-sided Student's t 
test. *P<0.05, “Р<0.01, ***P<0.001 vs. WT. Source data are available online. 


Definitions: NA — not applicable; NT, not tested. 
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Residue BRS3 Expression BA1 BA1 Emax MK-5046 pECso MK-5046 
Number mutants (%ofWT) pECsotS.E.M. +4S.E.M.(%WT) +S.E.M. EmaxtS.E.M.(%WT) 
à WT 1002060  -7.81+0.10 98.95+0.93 -9.340.08 99.6125.08 
253 196А 56.1440.57 -6.5340.06"* 194.35%2.509% -7.51%0.08%  169.5545.24** 

2.64 H107A 49.662225 -7.8220.08 7726+1579* MI MI 

3.29 5124А 88.2740.52 -7.980.03 121.41+172> -9.52+0.12 99.31+31.63 

3.29 5124р 46.9640.90  .7.833020 ვ.ვვანქ2“ -845:047 4727%11.52%% 

332 RI27A 7854149 -7.540.23 27.57+18 50“ -9.12+0.25 30.0944.74““ 

332  R127Q 88.5541.84  -8.0620.11  63.09%13.399% -9.06%016 53.85+10.57** 

333 1128А 69.7141.14 NA NA NA NA 

3.36 — V131A 49.8941.91 -6.9840.06* 211.09%33.36%% ^ -8.0440.15  220.25426.06** 

4.60 182A 63.4342.10 -7.13+006- 74.21Һ5,5%0 -4.93+3.56* NA 

ECL2 F200A 51624189 -6.98%0.79 110.6+38.69 MI MI 

ECL2 52054 61144037 -7.59:010  100.11+1874 -9.460.14 100.42+6.74 

ECL2 S205P 119.29+0.27 -7.3340.12  3057+41493**  -9.0640.34 —— 344748.97*** 

542 C221A 40.6240.58 -7.3840.18 129.13223.94 -867%021 96.99+4.55 

6.52 Н288А 34.0920.29 -7.1420.15** — 127.01418.88' -6.90.07*" NA 

6.55 Y291A 93.1340.58 -522%012%% NA -6.7730.21*** NA 

6.58 Н294А 94444120 -6.1540.16* 114.1413.68* NT NT 

6.58 H294R 76.8341.04 -8.3420.11* 110.127.27* NT NT 

ECL3  D303A 60114158 -7.3420.28 ^ 45.94+16.74"* NT NT 

7.31  H308A 82840177  -7.3120.27* 9707+1795 NT NT 

732 F309A 98.7442.23 -7.1740.06“ 109.52224.06 NT NT 

7.39  R316A 33.2540.30 -5.8840.04** NA -9.6240.04 176.84%25.63%" 

743 ЕЗ20А 43.7541.26 -6.874007“ — 45.3145.73** — -84940.558 | 40.2845.61““ 

E „ . 14.4140.08 -8.06+004 151.15413.19%* -9.65+0.11 89.3+3.87 
NMBR-WT  73.17+273 -8.6640.11** 285.63%35.44%% NA NA 
GRPR-WT 42.71+0.73 -8.8540.10"* 200.2%37.77%% NA NA 


1421 Table S5. Sequence alignment of the key residues in sodium site, DRY motif, 
1422 РУ(І)Е motif, toggle switch and NPxxY motif, as well as residues involved іп 


1423 disulfide bond formation in bombesin receptors. 


Sodium site Disulfide bond DRY motif PV(I)F motif toggle switch NPxxY motif 

Residue position 2.50 3.39 3.25 464 534 349 3.50 3.51 5.50 3.40 644 644 648 7.43 749 7.50 7.54 

NMBR D § C Т H D R Y P V E F уу Е N Р hé 

GRPR D S с 8 H D R Y Р V L Е VV F N ჩ né 

BRS3 D § C 8 L D R Y P V F F VV Е M P V 

Class A conserved D 8 C X x D R Y ჩ V(I) L F уу Е/Ү N P Y 
1424 
1425 
1426 
1427 
1428 
1429 
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